from django.contrib.gis.gdal import *
from myapp.models import GeoModel

ds = DataSource('/path/to/foo.shp')
lyr = ds[0]

for feat in lyr:
  lat = feat.get('LAT')
  lon = feat.get('LON')
  gm = GeoModel(point='POINT(%s %s)' % (lon, lat))
  gm.save()
